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One of the most fundamental quantities associated with 
polymer translocation through a nanopore is the translocation 
time t and its dependence on the chain length N. In a recent 
Letter, Wolterink et al. IUJ present new results for unbiased 
translocation using a 3D lattice model within the Monte Carlo 
(MC) simulation method. According to their scaling argument 
t ~ N 1+2v F(b/ R g ), where R g is the radius of gyration, b the 
pore width, the Flory exponent v = 0.75 in 2D and 0.588 in 
3D, and numerically the scaling function F(x) ~ a; ' _0 - 38±0 - 08 
for x — > 0. This leads to r ~ tv 2 - 40±0 08 , which is in con- 
tradiction with the established result by Chuang et al. |0] that 
r ~ N 2u+1 , which gives 2.50 in 2D and 2.18 in 3D. 

The original prediction that r ~ N 2v+1 [2] means that 
r scales with N in the same manner as the Rouse relax- 
ation time of the chain except for a larger prefactor, tr ~ 
R 2 g /D ~ N 1+2v , where the diffusion coefficient D ~ 1/7Y 
within Rouse dynamics. This result was recently corroborated 
by extensive numerical simulations based on the Fluctuating 
Bond (FB) [3] and Langevin Dynamics (LD) models with the 
bead-spring approach jj, 0], where r was found to scale as 

jy2.50i0.01 in m 

To resolve the apparent discrepancy, we have analyzed the 
approach of Wolterink et al. [1] and performed additional 
high-accuracy numerical simulations using the FB model with 
MC dynamics in 2D, and atomistic MD simulations us- 
ing the GROMACS [6] simulation engine in 2D and 3D. As 
in the MC and LD methods, explicit solvent hydrodynamics 
were excluded in our GROMACS simulations. We find that 
within numerical accuracy, both the FB and MD methods give 
the same scaling results in 2D, and thus here we present data 
as obtained using GROMACS within the bead-spring model. 
GROMACS is currently one of the most commonly used pro- 
grams in soft matter and biophysical simulations, and has also 
been used extensively by some of us in various problems (see 
e.g., Ref. 7 and references therein). 

Our main results are summarized in Fig. Q] where we find 
that r ~ 7V 2 - 44±0 03 in 2D and r ~ A 2 - 22±0 - 06 in 3D in 
complete agreement with Ref s . I^L 0> I^L ruling out results of 
Ref. [lj within the accuracy of the data. In particular, we find 
that the scaling function F(b/R g ) becomes independent of N 
for large N, as shown in the insert of Fig. 1 . Thus, there is no 
discernible correction to the exponent 1 + 2v. 

To verify our results independently, we also computed the 
squared change of the translocation coordinate s(t), where we 
observed sub-diffusive behavior ((As(t)) 2 ) ~ t a , with a — 
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FIG. 1: Scaling of escape times from GROMACS data. The slopes 
are 2.44 ± 0.03 (2D) and 2.22 ± 0.06 (3D). The insert shows the 
scaling function F(b/R g ) (~ t/N 1+2v ) as a function of N. 



0.807 ± 0.002 in 2D and a = 0.910 ± 0.002 in 3D. This again 
agrees with Chuang et al. [2] who obtained a = 2/(1 + 2u), 
which gives 0.8 in 2D and 0.92 in 3D. 

Further theoretical support for the exponent 1 + 2v comes 
from two other independent studies. In Ref. |3| it was ana- 
lytically predicted and also numerically confirmed that r ~ 
(R g + L) 2 / D for a pore of length L, resulting from the fact 
that the mass center of the polymer moves a distance of L . 
For long pore L ^> N we have r ~ NL 2 3> A 3 , which is 
longer than the reptation time of the chain ~ A 3 . 

In the second study a different approach was used: the start- 
ing point was driven translocation in which a pulling force F 
is acting on one end of the chain. In that case, r ~ N 2v+1 
can be analytically derived in the limit F — > using well- 
established scaling functions for polymers under tension 

To conclude, all the above independent results, analytical 
and numerical, and in particular the behavior of the scaling 
function F(b/R g ), confirm the result r ~ N 2w+1 , and invali- 
date the lattice model results of Ref. 1. The apparent discrep- 
ancy may be due to the artificial dynamics of the lattice model 
of Ref. 1 , which is based on the repton model for a single rep- 
tating polymer with the addition of sideways moves and repta- 
tion moves; although their lattice model works well for static 
properties in bulk solutions [9], the dynamics in the presence 
of a narrow hole may pose problems as the hole adds a new 
length scale to the problem. An off-lattice approach, such as 
used here is guaranteed to be free of any artifacts. 
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